Probing dark energy perturbations: 
the dark energy equation of state and speed of sound as measured by WMAP 
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We review tlie implications of having a non-trivial matter component in the universe and the 
potential for detecting such a component through the matter power spectrum and ISW effect. We 
adopt a phenomenological approach and consider the mysterious dark energy to be a cosmic fluid. 
It is thus fully characterized, up to linear order, by its equation of state and its speed of sound. 
Whereas the equation of state has been widely studied in the literature, less interest has been 
devoted to the speed of sound. Its observational consequences come predominantly from very large 
scale modes of dark matter perturbations {k < 0.01/iMpc~^). Since these modes have hardly been 
probed so far by large scale galaxy surveys, we investigate whether joint constraints that can be 
placed on those two quantities using the recent CMB fluctuations measurements by WMAP as well 
as the recently measured CMB large scale structure cross-correlation. We find only a tentative 
Isigma detection of the speed of sound, from CMB alone, < 0.04 at this low significance level. 
Furthermore, the current uncertainties in bias in the matter power spectrum preclude any constraints 
being placed using the cross correlation of CMB with the NVSS radio survey. We believe however 
that improvements in bias through improved survey scales and depths in the near future will improve 
hopes of detecting the speed of sound. 



I. INTRODUCTION 

With the recent unveiling of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) results, mea- 
suring the Cosmic Microwave Background (CMB) 
anisotropy 1], the on-going supernovae searches 3 
and the upcoming completion of the Sloan Digital Sky 
Survey, amongst others, we are seeing a wealth of pre- 
cision observational data being made available. To a 
great extent the standard A-CDM scenario fits the 
data well . 

However, the WMAP data might suggest that some 
modifications to the standard scenario are needed. 
One possible hint at required modifications is the 
deficit of large scale power in the temperature map, 
and in particular, the low CMB quadrupole whose 
posterior probability is less than a few hundredth (see 
e.g. for possible interpretations 0, Q and for a 
discussion of this number). One possibility is that 
this lack of large scale power might point to some 
particular properties of the dark energy. The domi- 
nant contribution to fluctuations on these scales is the 
Integrated Sachs Wolfe (ISW) effect which describes 
the fluctuations induced by the passage of CMB pho- 
tons through the time evolving gravitational potential 
associated to nearby {z < 5) large scale structures 
(LSS). One property we expect of dark energy is that 
it suppresses the gravitational collapse of matter at 
relatively recent times, which in turn suppresses the 
gravitational potential felt by the photons, thereby 
leaving a signature in the ISW correlations. Since 
this signature is created by the time evolving poten- 
tial associated with relatively close LSS, it should be 
intimately correlated with any tracer of LSS 0, 13 • 
A positive detection of such a cross correlation using 



WMAP data, assuming a cosmological constant as the 
dark energy, has recently been measured j2, ^| . 

However the underlying cause of the dark energy is 
still unknown; and such observational inferences offer 
rich prospects for guiding and leading the theoreti- 
cal effort. A wide variety of models have been pro- 
posed to explain observations, from the unperturbed 
cosmological constant, to a multitude of scalar field 
quintessence and exotic particle theories (see [l^ for 
a review). 

Much effort has been put into determining the equa- 
tion of state of dark energy, in an attempt to con- 
strain and direct theories. Since the equation of state 
affects both the background expansion and the evo- 
lution of matter perturbations there are a wealth of 
complementary observations available (again see |12| 
and references therein). 

An equally insightful, but less investigated, charac- 
teristic of dark energy is the speed of sound within it. 
This does not affect the background evolution but is 
fundamental in determining a dark energys clustering 
properties, through the Jeans scale. It will, therefore, 
have an effect on the evolution of fluctuations in the 
matter distribution. 

Following the papers laying the foundations for cos- 
mological perturbation theory [isLIT^ . the effect of the 
speed of sound on observables was considered in more 
detail: for the CMB and large scale structure 0, 
and in the context of lensing ■ Observational im- 
plications of the speed of sound in a variety of dark 
energy models have also recently been discussed: for 
example for k-essence , condensation of dark 

matter j23| and the Chaplygin gas, in terms of the 
matter power spectrum 21, 22, 2^ and combined full 
CMB and large scale structure measurements 0, . 
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Minimally coupled scalar field, quintessence mod- 
els, commonly have a non-adiabatic speed of sound 
close or equal to unit y ( in units of c, the speed of 
light), (see for example |2g). By contrast however, the 
adiabatic Chaplygin gas model, ( e.g. motivated by a 
rolling tachyon '27]) has a speed of sound directly pro- 
portional to the equation of state, both of which are 
approximately zero up until late times when the dark 
energy starts to dominate. It is conceivable therefore 
that distinctions between such models might be able 
to be made through the detection of a signature of 
the dark energy speed of sound: in the large scale 
ISW correlations, and in the cross correlation of the 
CMB with the distribution of large scale structure § . 

In section^we briefly review parameters describing 
a general fluid and the issues that arise when estab- 
lishing a fluid's speed of sound. In section HTll we de- 
scribe the implications of equation of state and speed 
of sound on perturbation evolution in the fluid and 
CDM. We consider a toy model with slowly varying 
equation of state and sound speed applicable in a wide 
variety of minimally coupled scalar field theories. In 
section ITVI we discuss the potential for detecting the 
speed of sound using late time perturbation evolu- 
tion, in the ISW effect, through the autocorrelation 
of the CMB temperature power spectrum. In section 
IVI we extend the analysis to the cross correlation of 
the WMAP CMB data with distribution. Finally, in 
section IVll we summarize our findings. 



II. THE SPEED OF SOUND WITHIN 
GENERAL MATTER 

For a perfect fluid the speed of sound purely arises 
from adiabatic perturbations in the pressure, p, and 
energy density p and the adiabatic speed of sound, 
is purely determined by the equation of state 
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where the subscript i denotes a general specie of mat- 
ter, where dots represent derivatives with respect to 
conformal time and where Ti. is the Hubble constant 
with respect to conformal time. 

In imperfect fluids, for example most scalar field 
or quintessence models, however, dissipative processes 
generate entropic perturbations in the fluid and this 
simple relation between background and the speed of 
sound breaks down and we have the more general re- 
lation 
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In order to establish the speed of sound in these cases 
we must look to the full action for the fluid described 



often through the form of an effective potential. In 
this case, the speed of sound can be written in terms 
of the contribution of the adiabatic component and 
an additional entropy perturbation and the density 
fluctuation in the given frame 5i [T^. 
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Fj is the intrinsic entropy perturbation of the mat- 
ter component, representing the displacement be- 
tween hypersurfaces of uniform pressure and uniform 
energy density. In this paper, we are solely interested 
in probing the intrinsic entropy of the dark energy 
component. It is worth noting that in a multi-fluid 
scenario, in addition to the intrinsic entropy pertur- 
bations denoted by F^, further contributions to the to- 
tal entropy perturbation of the system can arise from 
the relative evolution of two or more fluids with differ- 
ent adiabatic sound speeds , an d through non-minimal 
couphng (see for example |2^'). 

Whereas the adiabatic speed of sound, Cai, and F^ 
are scale independent, gauge invariant quantities, Csi 
can be neither. As such the general speed of sound 
is gauge and scale dependent and issues of preferred 
frame arise. Looking at equation Q, since the fluid i 
rest frame is the only frame in which Si is a gauge 
invariant quantity, this is the only frame in which 
a matter component's speed of sound is also gauge- 
invariant. 

A useful transformation fljj relates the gauge- 
invariant, rest frame density perturbation, Si, to the 
density and velocity perturbations in a random frame, 

Si and 9i, 



s,^s, + 3nii + w,)-^ 
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where we assume that the component is minimally 
coupled to other matter species and henceforth dark 
energy rest frame quantities are denoted using a cir- 
cumflex (') . 

Using equations |SJ) and lO we can rewrite the pres- 
sure perturbation in a general frame, Spi, in terms of 
the rest frame speed of sound. 



Spi = cliSpt + 3H{l + Wi){cl 



III. THE SPEED OF SOUND AND 
PERTURBATION EVOLUTION 

Herein, we use the synchronous gauge and follow 
the notation of CDM rest frame quantities are 

denoted Si and 9i, while dark energy rest frame quan- 
tities use the circumflex {"). The two are related by 
equation 
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FIG. 1: Relative suppression of dark energy perturbations 
to those in CDM as one increases ci from to 1 (top to 
bottom in both panels) for w = —0.8. The top panel shows 
perturbations in CDM rest frame, which due to entropy 
perturbations can become negative. The bottom panel 
shows perturbations in dark energy rest frame which are 
always positive but for w < are a fraction of the CDM 
perturbations, h and 0,^ are fixed so as to fit the WMAP 
constraints on Q.bh^ ,0.ch? and angular diameter distance 
to last scattering 



The energy density and velocity perturbation evo- 
lution of a general matter component in the CDM rest 
frame is given by 
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This set of equations illustrate clearly that linear per- 
turbations can be fully characterised by two numbers 
(and their potential time evolution): the equation of 
state and the rest frame speed of sound. 

Let us now consider a toy model with a general 
fluid in which the time variation in w and is small 
in comparison to the expansion rate of the universe so 



that we can model it with constant w (i.e. 



and 



Cg. Such models are not impractical and can be used 
as the basis for comparison with scalar field theories 
such as those with scaling potentials and Chaplygin 
gases during the radiation and matter dominated eras. 

In the matter dominated era, ignoring baryons for 
simplicity, CDM density perturbations are affected 
by the speed of sound of dark energy (denoted 'x') 



FIG. 2: Evolution the ISW source term, in comparison to 
the cS, — scenario, with w — —0.3, for = 0.25, 0.5, 0.75 
and 1 from bottom to top. h and 0,^ are fixed as in Fig. 
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through the relation, 
37^2 

Sc + Ti-Sc ^c^c 



(10) 



(1 + 3cl)S.^ + m{l + w){cl - w)- 



In the radiation and matter dominated eras the ex- 
pansion rate obeys an effective power law a oc r™ and 
we find that the evolution equations admit a solution 
of the form 6c oc S^/ {I + w) oc r^, and 0^ oc r^. 

Equation ((SJ shows how w and affect the relative 
size of dark energy and CDM perturbations. In Fig. 
n we see that, as one expects, increasing the speed 
of sound accentuates this suppression produced by re- 
ducing w, lowering Sx/Sc- Note that the dark energy 
density perturbation is well defined (remaining posi- 
tive) in the rest frame, while the transformation into 
the CDM rest frame can make dx negative; this is just 
a foible of the frame one chooses, however. 

The presence of dark energy perturbations leaves a 
w and dependent signature in the ISW source term. 
This can be written in terms of the time variation 
of the anisotropic stress and the rest frame density 
perturbations of each matter component. 
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Since Sx is suppressed in comparison to Sc, the dom- 
inant contribution to the ISW will come from the 
CDM perturbations. Subsequently it will be suppres- 
sion of these (in comparison to a Qc — ^ scenario) 
through the effect of the dark energy speed of sound 



12 3 

log 1 



-0.8 -0.6 -0.4 -0.2 
w 



FIG. 3; CMB TT spectra for w = -0.3 (top panel) and 
w = —0.9 (bottom panel) with cl — and 1, all other 
parameters fixed to give the best fit at smaller scales. All 
spectra are normalized to Cs7 for comparison. 



FIG. 4: Likelihood contour plot for the dark energy com- 
ponent in the w — Cs plane showing 1, 2 and 3a contours 
(heaviest to faintest lines respectively) after marginalising 
over the power spectrum normalisation . 



and equation of state that will leave a signature in the 
ISW. Fig 2. shows how as one increases the ISW 
effect increases. 



IV. CONSTRAINTS USING WMAP 
TEMPERATURE FLUCTUATIONS AT LARGE 
SCALES 

In this section, we investigate the joint constraints 
on the equation of state and the speed of sound that 
can be infered using the CMB temperature power 
spectrum. As was discussed earlier, the main effect 
of a sound speed smaller than the speed of light will 
be felt at late times and large scales and will thus 
only affect the very large scale CMB temperature fluc- 
tuations arising from the late ISW elfect, for which 
WMAP already provide us will full sky cosmic vari- 
ance limited measurement 1]. The Fourier compo- 
nent of the fluctuations arising from the ISW effect 
are given by (for this equation only, we ignore the well 
known issues related to the sphericity of the observed 
sky). 
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where is the square of the speed of light, Hq is 
the Hubble constant today, (5°(fc) — Sc{z — 0,fc) 
and is the fractional energy density in matter 
(CDM-|-baryons) today and g{z, k) = {1 + z)D{z, k) 



where D{z,k) is the linear growth factor given by 
5c{z^ k) = D{z, k)6'^{k). In the limit that w tends to -1 
D(z, k) is scale independent and can be approximated 

by 
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however this approximation does not provide the de- 
gree of precision that is required for w > —1. even 
in the absence of dark energy perturbations 6u\. Be- 
cause of this and in order to factor in the late time 
effect of dark energy perturbations and their scale de- 
pendence, for 7^ 0, we explicitly calculate the linear 
growth function, D{z, fc; c^, w) for each model. 

Note that the effect of the speed of sound comes in 
solely through the value of (5° while the equation of 
state affects both 5° and the linear growth factor. 

The associated auto power spectrum is given sub- 
sequently by 
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dk 



^(fc)^(fc'))j|(fc'?)-(16) 



In order to probe solely the effect of on Cf^ , 
we will compare to WMAP observations a family of 
models lying along the angular diameter degeneracy 
surface present in the CMB spectrum. To do so, we 
keep ^Ich-^ = 0.135 and Ob/i^ — 0.0224 to be consistent 
with the WMAP best fit '3] and choose w and h such 
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that the angular diameter distance to last scattering 
is the same. Other parameters correspond to the best 
fit model of Q (table 7). Doing so, we ensure that 
only the large scale correlations vary with each model 
and that in all other respects they fit the WMAP data 
well. Note that we also have vary slightly the overall 
amplitude due to the change in the first peak height 
ISW plateau ratio. We consider w values between 
and -1 and between and 1. Given this grid of 
model, we can then deduce easily the likelihood of the 
data using the publically available code provided by 
the WMAP team (sj , from which we can deduce some 
joint constraints on w and c^. 

In Fig. 13 we show the variation of the CMB TT 
power spectrum as one varies from to 1 for a 
model with w = —0.3 and -0.9. One can see that in- 
creasing increases the suppression of the CDM per- 
turbations and therefore increases the power on large 
scales. The effect decreases though as one decreases 
w; at low w, the suppression due to the equation of 
state itself will generate a dominant ISW effect on top 
of which a subdominant contribution from is then 
superimposed. Those results agree with the one ob- 
tained in [TtI I. 

In Fig. 21 we show the likelihood plot from the 
WMAP data in the w — c^ plane. The low quadrupole, 
and other low i Cg's lead to a value of < 0.04 being 
preferred by the data, at the la level, although as one 
moves to lower w the ability to distinguish between 
different values of the speed of sound disappears, be- 
cause of cosmic variance. 

The cosmic variance thus limit our ability to con- 
straint the dark energy speed of sound using temper- 
ature Ci only. However, given the fact that all the 
constraints comes from the ISW effect, it is natural to 
consider the cross correlation the CMB with the large 
scale distribution of matter near us, correlation that 
is a direct probe of the late ISW. In theory then, this 
might give us a better and different probe into c^, so 
that both should be combine eventually. We consider 
this in the next section. 



V. CONSTRAINTS USING CMB AND 
LARGE SCALE STRUCTURE CROSS 
CORRELATION 

As stated earlier, the dark energy affects very large 
scale modes of dark matter density perturbations. As 
shown on Fig. [S] , those modes are outside the range 
of current wide field galaxy surveys. For example, the 
SDSS measured the galaxy power spectrum down to 
~ 0.01 /iMpc^i "only" (see e.g. [H). Full sky survey 
exists though, but their particular properties and in- 
trinsic limitations restricted their use as direct probes 
of the matter power spectrum at those scales. For 
example, the NRAO VLA Sky Survey (NVSS) [H 
encompass such a wide variety of objects that the dif- 
ficulties in modeling the biases at stake prevented its 
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FIG. 5: Matter power spectra spectra COBE normalised 
Top panel: w = -0.9 (dashed) and -0.3 (full) with = 
and 1 (top and bottom lines respectively ) 
Bottom panel: = 0, it; = —1, —0.75, —0.5, —0.25 from 
top to bottom (at large scales on both plots) 



usage to directly measure (dark) matter density fluc- 
tuations at those scales and infer this way any precise 
cosmological constraint. However, their use in con- 
junction with large scale CMB fluctuation measure- 
ments allows us to circumvent somehow this difficulty. 

Indeed, within a given cosmological model, one can 
look at the surveyed objects as a simple linearly bias 
tracer of dark matter perturbations, a reasonnable ap- 
proximation on those very large scales. By measuring 
the auto-correlation function (ACF) of those objects 
on those scales, one can infer the model dependant ef- 
fective bias for this composite population. Since this 
population trace the large scale gravitational poten- 
tial, it should correlate with the CMB fluctuations 
induced by the same potential through the ISW ef- 
fect. The angular dependance of this cross-correlation 
function (CCF) and its amplitude both depend on the 
tracer properties (bias and redshift distribution) and 
on the particular cosmological model considered. In 
particular, we would expect an important dependacy 
on the dark energy properties which drive the evo- 
lution of the universe at those late times and large 
scales. Modeling the tracer properties, we can thus 
in principle constrain the cosmology. This has been 
advocated flrst in [^, studied in details in 0, 
and performed effectively using as a tracer the NVSS 
sources IlllSlal, HE AO- IX- ray sources [s^, HI] or 
APM galaxies 

So far, this correlation has been probed to prove 
the very existence of dark energy and to constrain its 
overall density. We here extend this approach and try 
to investigate the potential constrains on its very per- 
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FIG. 6: Predicted ACF for the NVSS radio source catalog 
as a function of angular separation, 9,, in degree. The lin- 
ear bias factor has been accounted for. The various mod- 
els correspond, in orange, to w = —0.9 and Cg = 0.9,0.1 
(respectively in triple dot-dashed or dot-dashed lines), or, 
in blue, to w = —0.3 with — 0.9,0.1 (respectively in 
triple dot-dashed or dot-dashed lines). The black line cor- 
responds to our fiducial ACDM model. Measurements in 
red are from ^10^1 . 
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FIG. 7: Predicted CCF for the NVSS radio source catalog 
as a function of angular separation, 6, in degree. The lin- 
ear bias factor has been accounted for. The various mod- 
els correspond, in orange, to w = —0.9 and = 0.9,0.1 
(respectively in triple dot-dashed, dot-dashed lines), or, 
in blue, to w = —0.3 with — 0.9,0.1 (respectively in 
triple dot-dashed or dot-dashed lines). They largely over- 
lap. The black line corresponds to our fiducial ACDM 
model. Data points from llQjl are in red as well as (corre- 
lated) error bars. 



turbative properties, i.e. jointUy its equation of state 
and its sound speed. We will use as a data-set the 
ACF and CCF measurements of performed using 
the NVSS catalog and the WMAP l-year maps. For 
the sake of simplicity we will follow the same notations 
that we recall brievly. 

The measurements of fluctuations in the nearby 
matter distribution from measuring the radio sources 
distribution can be expressed in terms of the fractional 
source count perturbation given by, 



(17) 



where hr is the linear bias in the matter distribu- 
tion, Nq is the mean source count per pixel (147.9 
for l.Sdeg square pixels used in ^3)' ^'^'^ dN/dz is 
the normalised redshift distribution of galaxies, such 
that [ (dN I dz)dz = 1. For the latter we adopt the 
model of |33 |. 

The dimensionless two point correlation function 
between two quantities X and Y with background val- 
ues Xq and Yq in positions h and fi' in the sky is given 

by 



5X ,^JY 
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= ^y\2l + l)Cf^Peico.s9) (18) 

For the fractional source count and CMB temperature 
cross and auto correlations 
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(19) 



^A2(fc)/f (fc)/f (fc) , (20) 



A^{k) = {S",{k)5°{k')) = S{k - k')k^P{k)/2n^ (21) 

and where the filter functions for the temperature and 
number count fluctuations respectively are given by 



_ r dg{z,k)^ 

c^k"^ jr. dz 



f^{k) = h 



j,{kri{z))dz (22) 



dN 
dz 



D{z,k)j,{krj{z))dz , (23) 



where D{z, k) has been previously deflned in section 
M 

The results of this approach are summarised in fig- 
ures El and [3 As was used in section Hvl we consider 
here a family of models lying along the CMB angu- 
lar diameter degeneracy surface, and move along it by 
varying simultaneoulsy w and h. For each of these 
background, we then consider various and compute 
the linear transfer function using a modified version 
of the CMBfast gO;] or CAMB ^] softwares . For 
a sample of those models, we plotted both the pre- 
dicted ACF, C^'^{9), from which we infer the bias 
(figureini) and the subsequent predicted CCF, C^'^(9) 
which one can compare with the measurements of |lO| 
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(figure [TJ. Obviously, even if this signal does indeed 
have some dependance with regards to the dark en- 
ergy pertubations parameters, w and c^, most of the 
effect is absorbed in the bias measurement illustrated 
in figure The uncertainty in the bias is the main 
hindrance in placing constraints on both the equation 
of state and speed of sound, so that current data do 
not allow this correlation is not measured well enough 
to constrain our models. We obtain a bias ran ge a nd a 
bias evolution similar as the one obtained by |1CI |. i.e. 
br ~ 1.4 for our fiducial ACDM model, and br tends 
to decrease when w increases (in range comprised be- 
tween 1.4 and 2.2). Note that the plotted error bars 
are heavily correlated. The knowledge of the corre- 
lation matrix computed in [Tol | allows us however to 
compute a well defined statistic. Note also that 
those error bars include cosmic variance only but do 
not take into account the shot noise in the radio source 
catalog. It would be important to consider the uncer- 
tainties in the bias measurement that propagates di- 
rectly in the signal prediction but we leave this issue 
for future work. Nevertheless, given the fact that all 
our models lie within the one sigma error bars, we will 
not illustrate here by a contour plot those weak joint 
constraints on w and c^. 

A natural and important question that arises at this 
level concerns the future prospect for the measure- 
ments of this correlation, depending on the used LSS 
tracers as well as the survey considered. Although 
some studies have already been done 0,13 H^, more 
specific investigations are necessary. In particular an 
independent measurement of the bias, along with im- 
proved scale and depth of survey will all contribute to 
vastly improving on the current observational uncer- 
tainties. 



VI. CONCLUSIONS 

We have reviewed the effect of the speed of sound 
of dark energy on CDM and dark energy perturba- 
tions. While a positive dark energy speed of sound 
suppresses the CDM perturbations, it is the deviation 



from adiabaticity, in combination with the equation of 
state that determine the degree of suppression of the 
amplitude dark energy perturbation in comparison to 
those of the CDM. 

We have found the CMB large scale temperature 
fluctuations, dominated by the ISW effect, are a 
promising tool to measure the speed of sound. The 
suppression of CDM matter perturbations drive the 
late time ISW effect. 

From the auto correlation of the WMAP data with 
itself we obtain a la "constraint" on the speed of 
sound < 0.04, using scenarios that minimise con- 
tributions to the likelihood on small scales (from the 
peaks) as much as possible by using well known degen- 
eracies to follow the WMAP best fit model as closely 
as possible. The main limitation in obtaining con- 
straints from the auto correlation is the cosmic vari- 
ance. 

We have also investigated the cross correlation of 
the large scale CMB with fluctuations in the nearby 
mass distribution using the NVSS radio source cat- 
alogue. We here again find that cosmic variance is 
a strong limitation and prevent us from placing any 
strong constraint in the w — cj. plane. 

However, since the potential of such an analysis 
might be unique in unveiling the mysteries of the 
dark energy, it is important to explore further out 
the prospect of future potential large scale probe of 
the gravitational portential and so of the ISW (LSST, 
PLANCK, CMBPOL). We have presented some esti- 
mates of prospective constraints that one might obtain 
from cross correlation of large scale probes with CMB 
however we leave this exploration for future work. 
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